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ABSTRACT 



Aims. We search for frozen water and its processing around young stellar objects (YSOs of class I/Il). We try to detect potential, 
regional differences in water ice evolution within YSOs, which is relevant to understanding the chemical structure of the progenitors 
of protoplanetary systems and the evolution of solid materials. Water plays an important role as a reaction bed for rich chemistry and 
is an indispensable requirement for life as known on Earth. 

Methods. We present our analysis of NAOS-CONICA/VLT spectroscopy of water ice at 3 /im for the T Tauri star YLW 16 A in the 
p Ophiuchi molecular cloud. We obtained spectra for different regions of the circumstellar environment. The observed absorption 
profiles are deconvolved with the mass extinction profiles of amorphous and crystallized ice measured in laboratory. We take into 
account both absorption and scattering by ice grains. 

Results. Water ice in YLW16A is detected with optical depths of between t = 1.8 and t = 2.5. The profiles that are measured can be 
fitted predominantly by the extinction profiles of small grains (0. 1 yum - 0.3 ;um) with a small contribution from large grains (< 10%). 
However, an unambiguous trace of grain growth cannot be found. We detected crystallized water ice spectra that have their origin in 
difi^erent regions of the circumstellar environment of the TTauri star YLW 16 A. The crystallinity increases in the upper layers of the 
circumstellar disk, while only amorphous grains exist in the bipolar envelope. As in studies of silicate grains in TTauri objects, the 
higher crystallinity in the upper layers of the outer disk regions implies that water ice crystallizes and remains crystallized close to the 
disk atmosphere where water ice is shielded against hard irradiation. 

Key words. Infrared: stars - Accretion disks - Astrochemistry 
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1. Introduction 

The importance of water (ice) to protoplanetary systems is jus- 
tified by the following properties: liquid water is a proper sol- 
vent for many minerals and organic molecules and allows the 
formation of complex organic molecules. Water is the reac- 
tion bed for the photochemical synthesis of both aromatic com- 
pounds (Bernstein et al. '2002^ and amino acids (Munoz Caro et 
al.|2002), even in its solid state. Therefore, rich deposits of wa- 
ter (ice) in protoplanetary disks are assumed to be a requirement 
for the formation of life as known on Earth. The importance of 
water ice is in addition emphasized by the assumption that ice 
planets, such as Neptune, form outside the snowline where wa- 
ter condenses on dust grains. The mass fraction of solid matter 
in protoplanetary disks abruptly increases at the snowline at least 
by one order of magnitude (Stevenson & Lunine[T988]l. 

The water ice band at 3.08 ;um, which is caused by a vibra- 
tional excitation of the OH bond in the water molecule, was 
primarily discovered in the envelopes of numerous deeply em- 
bedded protostars (Gillett & Forrest [19731 Gillett et al. [19751 
Merrill et al. jl976 ). Frozen water that can be found in dust par- 
ticles is the most frequent ice molecule in YSOs (Ehrenfreund et 
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al. '1996'01 The huge water ice deposit results from its molec- 
ular property: in contrast to nonpolar, volatile ice molecules, 
such as CO, CO2, N2, and O2 ice, H2O ice has a high subli- 
mation temperature of ~120K depending on the ambient pres- 
sure (Davis I2OO7I 1. The polarity of the water molecule is re- 
sponsible for its stable adsorption by condensation seeds such 
as silicate dust grains. Cosmic irradiation and chemical reac- 
tions could force the evaporation of water ice (Leger et al. 1983J 
Hartquist & Wolf 1990). However, these reactions are impeded 
by regions with strong shielding and low temperatures (Tielens 
& Hagen 1982 ). This is why the water ice band is commonly ob- 
served in the spectra of deeply embedded protostars (Boogert et 
al. 120041) . However, it is assumed that larger deposits of water ice 
exist in the optically thick regions of circumstellar disks (Fig. 1 
in Chiang et al. l2001l l. This assumption may be supported by the 
detection of the water ice band at 46.5 //m (Creech-Eakman et 
al. 2002) and 3.08 yum (Leinert et al. l200Tl Terada et al. lIOOTl ) in 
the spectra of several TTauri objects. 

The water ice band at 3.08 /vm in L bancQ is assumed to be 
an excellent indicator of the evolutionary status of circumstel- 



' Fourteen different ice compounds, such as CO, CO2, or CH4 ice 
are known. Ice molecules are weakly bounded just by van-der-Waals 
forces or by hydrogen bonds. A chemical compound that can be found 
in metals, minerals, and diamonds does not form. 

^ The L band covers a wavelength interval between 2.8 //m and 
4.2 /jm in the near-infrared range. 
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lar disks (van de Bult ll985T l. In a protostellar molecular cloud, 
water ice has a primarily amorphous structure, as many previ- 
ous studies have shown (e.g., Thi et al. 2006; Pontoppidan et 
al.|2005). As soon as the temperature exceeds a level of ~80K, 
water ice begins to crystallize (Hagen et al. |1981aD . This con- 
version can be detected by monitoring the profile of the optical 
depth, i.e., the logarithm of the absorption band: as the crys- 
tallinity increases, the profile narrows and its maximum shifts to 
longer wavelengths (e.g.. Smith et al. 119891 ). But the crystalline 
structure can also be destroyed, e.g., by hard irradiation (Kouchi 
& Kuroda [19901 ). In addition to crystallization, grain growth of 
ice particles can be inferred from a modification of the optical 
depth profile. An increase in the dust grain radius results in a 
higher scattering efficiency, a broadening of the absorption band 
at long wavelengths, and a shift of the maximum to longer wave- 
lengths. These changes are comparable to those related to the 
prominent 10/vm silicate band (e.g., Bouwman et al. 2001), but 
the origin of the spectral bands is different. While the emission 
band of silicate has its origin in the warm disk atmosphere, the 
absorption band of water ice originates in more embedded, cold 
regions closer to the disk midplane and far from the central star. 

The water ice band interferes with absorption bands of ad- 
ditional materials that are dissolved in water ice or form a 
complex compound with the water molecule: ammonia hy- 
drate H2O.N H3 (m inimum at 3.47 ;um: Mem ll et al. 119761 
Mukai et al. 119861 Dartois & d'Hendeco urt 120011 ). ammo- 
nia NH3 (2. 96 jum: d'Hendecourt et al. 119851 Dartois & 
d'Hendecourt 119851 ). methanol CH3OH (3.54//m: Dartois et 
al. [2001), hydrocarbons with the vibrational excitations of the 
CH, CH2, and CH3 bond (3.38 /jm, 3. 42//m, and 3.48 jum: Duley 
& WilHams 1984; Chiar et al. 120021) . and other chemical com- 
pounds (Ehrenfreund et al. 1996). The absorption bands of these 
molecules are narrow and interfere predominantly with the wa- 
ter ice band at longer wavelengths around 3.6 //m. They can 
therefore be distinguished from the water ice band. To study the 
chemical evaluation of ices in solar-mass systems, a survey in 
the infrared wavelength range where 41 low-luminosity YSOs 
were included was conducted by using the infrared spectrograph 
on the Spitzer Space Telescope (Boogert et al. 2008 ). 

The scientific scope of this study is to look for spatial vari- 
ations in the water ice band profile at 3.08 //m within a TTauri 
star to detect possible, regional differences in water ice evolu- 
tion. In Sect. |2] and [3] we scrutinize in detail our spectroscopic 
observations and data reduction, respectively. The derived water 
ice profiles are presented in Sect. |4] The ice profiles that could 
be derived are analyzed in Sect. I5.2l bv considering both the ab- 
sorption profiles of ice grains of different size and crystallization 
degree that are measured in laboratory. This well-established ap- 
proach of deconvolution is described in Sect. 15.11 The question 
of whether evolved, i.e., large and/or crystyllized ice grains may 
finally be found is discussed in Sect. I5.3l and l5.4l We conclude 
this study with a summary in Sect. |6] 

2. Observations 

Our observations with the near-infrared (NIR) spectrometer 
CONICA (Lenzen et al. 2003) were performed at the 8.2 m 
YEPUN telescope of the Very Large Telescope (VLT) in spring 
and summer 2006. We used the L54_1_SL mode where a 
wavelength range of between 2.60 yum and 4.20 /im is covered. 
Observations with CONICA are supported by the adaptive op- 
tics NAOS (Rousset et al. 120031 ). A linear, spectral dispersion of 
3.16 nm /pixel could be reached. The spatial resolution that could 
theoretically be reached was 0.054". Depending on the weather 



conditions during our observations, the pixel scale was ~0.12". 
The slit of the spectrometer had a width of 0.172" and a length 
of 28". 

We selected YSOs with quasi edge-on disks to be able to ob- 
serve optically thick regions that shield ice from hard irradiation. 
Images, polarimetric maps, and previous object models were 
used as selection criteria. A visual extinction of Ay > 3 mag 
is sufficient for an effective shielding of water ice against irradi- 
ation (Murakawa et al. 120001 ). 

A lunar occultation observation of YLW16A at 2yum 
(Simon et al. 1987) showed that the object is a single extended 
source and strongly inclined. Two brightness peaks are separated 
by -0.5" in a K band image obtaine d by N ICMOS at the Hubble 
Space Telescope (HST; Allen et al. 120021 ). These peaks are as- 
signed to the optically thin bipolar envelope. 

Table [Tlprovides an overview of our observations of the tar- 
get YLW16A and the corresponding standard star HR6070. 
The brightness of the standard star, an A V star, is known from 
the standard star catalog of van Bliek et al. ( 11996 ). For the prepa- 
ration of the observations, the L band brightness of the target 
was estimated in this study by considering the H band magni- 
tudes and the (H-K) and (K-N) colors of additional YSOs of the 
same star-forming region (Allen et al. 2002). 

Our target was observed using the nodding mode with a nod- 
ding angle of 10". Thus, the target was successively observed 
in different detector areas where the sky background could be 
eliminated by subtracting subsequent exposures. The chopping 
mode, which provides a more effective elimination of short-term 
variations in the sky, was not available in the observation mode 
used. 

The target was observed twice, with two orthogonal orien- 
tations of the slit. The observation sequence started with a slit 
orientation parallel to the rotational axis of the object (position: 
p-parallel). In a second exposure, the slit was rotated by 90° (po- 
sition: o-orthogonal). The photometric center served as a rota- 
tion center This procedure is assumed to allow the determination 
of the spatial distribution of water ice in the different regions of 
the circumstellar environment. 

In the appendix of this publication, we list additional obser- 
vations of YSOs within this program. These YSOs also belong 
to the p Ophiuchi region. 



3. Data reduction 

When reducing the spectroscopic data obtained with NAOS- 
CONICA, particular problems have to be taken into account. The 
specific steps of the reduction are therefore considered in detail: 

i. Bad pixels are localized and replaced by the median of the 
pixels within the surrounding 3x3 detector segment. 

ii. The background of detector, telescope, and sky are elimi- 
nated by the subtraction of successive exposures for different 
nodding positions. Owing to unavailability of the chopping 
mode that can clear short-term variations in the sky, telluric 
features cannot be completely removed in the object spectra. 
The resulting images are flatfielded. 

iii. The spatial and spectral dimensions of a spectroscopic ex- 
posure are assigned to the rows and columns of the image 
matrix, respectively, to a first approximation. A Gaussian 
function is fitted to each row of the array to determine the 
maximum of the spatial intensity distribution. The counting 
rates of the pixels within 4 x FWHM (full width at half max- 
imum) of each Gaussian function are added. Alternatively, 
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Table 1. Overview of our observations with NAOS-CONICA. Coordinates (Cutri et al. 120031) . L band magnitude, spatial FWHM, 
date, observing time (UT), as well as the mean airmass (AM) during the observations are listed. The object YLW 16 A was observed 
twice using two different slit orientations: parallel (p) and orthogonal (o) to the rotational axis of the YSO. The quantity Ttot is the 
total exposure time. Directly before or after the observation, the standard star HR 6070 was observed. 



object 



Q-(J2000) ^(J2000) L[mag] FWHM ["] 



date 



UT 



AM 



April 28'\ 2006 
April 28"', 2006 

August 10"', 2006 
August 10"', 2006 



6070 


161817.9 


-28 36 50 


4.807 


W16A(p) 


16 27 28.0 


-24 39 34 


6.7 


W16A(o) 


16 27 28.0 


-24 39 34 


6.7 


6070 


1618 17.9 


-28 36 50 


4.807 



0.5 
0.5 



00- 
26- 

45- 
51- 



02 
08 

00 
59 



1.0 
1.0 

1.1 
1.3 



2.0 s 
60 X 3 X 4.0 s 

26 X 50 X 2.8 s 
2.0 s 



the Gaussian functions are divided into seven single inter- 
vals resulting in spatially adjacent spectra. For a pixel scale 
of 0.12" and the source distance of 160pc, the single inter- 
vals have widths of about 20 AU, 20 AU, 30 AU, and 40 AU, 
respectively, symmetrically arranged of increasing width at 
increasing distance from the brightness maximum. The inter- 
val of the central region has a width of 20 AU, the intervals 
of the outermost regions have widths of 40 AU. 

iv. The wavelength calibration of the L band spectra acquired 
with NAOS-CONICA is difficult because observations of 
specific calibration lamps could not be performed in this 
observing mode. Intrinsic stellar emission lines such as the 
Pfy-Hne at 3.74;um (Wallace & Hinkle 2002) do not clearly 
emerge from the background noise. Therefore, the measured 
spectra were cross-correlated with the telluric features of a 
sky spectrum obtained from NSF/NOAOQ As the shapes of 
the asymmetrically curved telluric lines are not known, the 
accuracy of the wavelength calibration cannot be better than 
the accuracy provided by a pixel width. 

V. To perform subsequent flux calibration, the standard star 
HR6070 was observed in the immediate vicinity (~1°), di- 
rectly before or after the observation of the target. This en- 
sured that instrumental and the atmospheric transmission 
could be determined, simultaneously. The division by a 
black-body function By{T^f{) with T^ff = 9980 K for HR 6070 
(Allen 1200 0) helps to reduce the influence of the continuum 
of the standard star Considering the template spectrum of an 
FO V star in L band (Rayner et al. 120091 ). the template spec- 
trum of an AO V star up to 3.0;um (Pickles T998) and the 
amplitude of the absorption depth of the water ice feature, 
potential stellar lines are assumed to be neglected, i.e., we 
do not correct for potential absorption lines in the spectrum 
of HR6070. Figure [1] presents the normalized transmission 
of sky and instrument that were obtained from the observa- 
tion of the standard star HR 6070. Spectral shifts of the stan- 
dard star spectra on a sub-pixel scale allowed an improved 
elimination of the telluric lines in the spectrum of the target. 
To remove atmospheric residuals and data points affected by 
poor atmospheric transmission, sections of the spectra with 
less than 35% transmission are not considered in the anal- 
ysis. These vary from exposure to exposure depending on 
weather conditions. 

vi. By considering the photometric L' band spectrcQ and the fil- 
ter curves that were formerly used to perform the photomet- 
ric measurements of the standard star (van Bliek et al. 1996 ), 
a photometric flux calibration of the observed spectra of the 




^ The sky spectrum is available at 
|http : //www .esc. org/sci/f acilities/paranal/instrumentsT] 
'' The L' band covers the NIR range between 3.4 yum and 4.2yum. 



Fig. 1. Transmission of atmosphere and instrument derived from 
the spectrum of the standard star HR 6070. 

target is possible. We note that the resulting scale factor of 
the photometric calibration can vary by more than 15% dur- 
ing a night (Przygodda 2004) . 



4. Results 

Figure |2] shows the resulting L band spectra of the target for 
which the counting rates of 4xFWHM of the spectrum were de- 
termined (Sect. [3]). To provide spatial information, Fig.[3]shows 
the spectra derived from seven spatially adjacent stripes of the 
total spectra measured in orthogonal orientations. Spectra 1 to 
7 represent spectra derived from source regions with projected 
widths of 40 AU, 30 AU, 20 AU, 20 AU, 20 AU, 30 AU, and 
40 AU, respectively. 

The NICMOS/HST-images of the class-I-object YLW 16 A 
exhibit a bipolar envelope above and below its circumstellar disk 
(Allen et al. .20021 . Both bipolar components yield two separate 
spectra I and II in the (p)-orientation of the slit. 

The photometric flux can be compared to photometric mea- 
surements of the Infrared Array Camera (IRAC) onboard of the 
Spitzer satellite. We perform photometry on the pipeline mosaic 
available at the Spitzer archive^ The photometric calibration of 
IRAC is based on an aperture whose radius is set to be lOpixels. 
The flux amounts to 0.79 Jy at 3.6 /im after background sub- 
traction. In contrast to NAOS-CONICA, IRAC does not allow 
us to spatially resolve the geometrical structure of YLW 16 A. 
Therefore, the different spatial resolution powers of IRAC and 



pittp : //ir sa ■ ipac . caltech . edu/Mi ssions/spitzer . html | 
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NAOS-CONICA explain the differences between the L band 
fluxes. 

All spectra derived from our NAOS-CONICA observations 
display the broad water-ice band absorption feature with its min- 
imum at ~3.1fim. As the atmospheric transmission is low at 
wavelengths of <3.0//m (Fig. [T]), the water ice feature is trun- 
cated at this lower wavelength end of the spectra. 



5. Modeling of the water ice band 

5.1. Modeling approach 

By definition, the depth of the absorption band is the natural 
logarithm of the ratio of the measured flux Fobs to the continuum 
flux Fcont above the absorption band. It depends on the path L 
through the absorbing material, its (mean) extinction Qext, and 
the particle density p. Whittet (120021) derived the optical depth r 
of spherical particles with a mean radius a and a mean material 
density rj to be 



In 



^obs 



QextTTfl 



^ext''" V^ 

= Lp—- r- ^ LpX A-ext = ^0 > , [Ki X AText;,) ■ 



Al'ina^rj 



(1) 



The quantity /Cext represents the (total) mass extinction. The 
extinction results from a linear superposition of mass extinc- 
tions /Cext;/ of material / that are weighted with the coefficient Ki 
(Martin [T9781 Solc ll980l l. The quantity Kq = Lp is a constant 
for each object. 

To determine the ice components that eff'ectively contribute 
to the 3 /im-absorption band, extinction profiles of ice grains of 
different size and crystallinity are linearly combined to fit the ice 
profile. As a first approximation, scattering effects can only be 
neglected if the condition 



2na 1 
A \m\ 



(2) 



is fulfilled (Dartois et al. I2001I I. where the quantity \m\ is the 
complex refractive index. For \m\ w 1 .35 (Querry et al. II969I I and 
A - 3.5 ;um, we obtain a < 0.41 yum. This result justifies the con- 
sideration of scattered radiation. Furthermore, the bipolar enve- 
lope of the object YLW 16 A appears in scattered light in the cor- 
responding K band images obtained with NICMOS/HST. In con- 
trast to the absorption profiles, the extinction profiles are asym- 
metric, i.e., more bell-shaped at longer wavelengths (Dartois & 
d'Hendecourtl200Tll. 

For our fitting routine, we assume opacities of entirely amor- 
phous {T - 10 K) and crystallized {T > 140 K) water ice. The 
ice grain radii a are 0. 1 //m, 0.3 //m, 0.5 fim, 0.8 //m, and 1 .5 //m. 
In this study, grains with a ^ 0.5 //m and a < 0.5 ;um are called 
large and small grains, respectively. A restriction to two grain 
radii (e.g., 0.1 /vm and 0.8 //m) only increases the deviation be- 
tween model and measurement but does not modify the result- 
ing mass ratios Kj. Figure |4] shows the mass extinction Ke^t(A) 
of amorphous and crystallized water ice used in our fitting ap- 
proach (Eq. [U. Assuming spherical, compact ice particles, the 
mass extinction curves are derived using the program MIEX 
(Wolf & Voshchinnikov |2004d . which is based on Mie scatter- 
ing. This program calculates the extinction profiles by consid- 
ering the grain size and refractive indices (Schmitt et al. II998I 
Dartois & d'Hendecourt l2001l l. 



An unambigous determination of the continuum flux F^om 
above the absorption band is difficult. Owing to reduced at- 
mospheric transmission (Fig. [TJ, parts of the spectra at short- 
wavelengths are missing. The determination of the continuum 
at longer wavelengths (>3.6yum) is also difficult because some 
telluric lines could not be removed entirely during data reduc- 
tion, producing noisy spectra. In previous investigations of the 
water ice band (Thi et al. 120061 Dartois et al. 2002 ; Boogert et 
al. I2008I I. a Planck function By(T) or a spline function was fit- 
ted to the continuum and - if available - to photometric data 
points in adjacent spectral bands to determine the continuum. In 
addition, studies were performed where the spectral energy dis- 
tribution of the entire infrared wavelength range was fitted using 
disk models (Pontoppidan et al. 120041 ). Assuming a spatial reso- 
lution of 0. 12" (-20 AU) for our NAOS-CONICA observations, 
which is smaller than the true spatial extension of the object, 
only a fraction of the total flux of the disk is detected. Therefore, 
when comparing with IRAS observations (Sect. |4|, photometric 
measurements do not have to represent the absolute level of the 
L band spectra measured with NAOS-CONICA. 

Apart from either the Planck function By(T) or a spline 
function, we use a straight line to determine the continuum. 
The straight line is rotated successively around different spectral 
points resulting in different optical depths (Eq.[T]). The spectral 
points about which the line is rotated are the sampling points 
at 3.6 /im and 3.8;um, respectively. The latter point is used be- 
cause the water ice band and other compounds (Sect. [TJ do not 
contribute to the spectrum around 3.8 yum. The angular step used 
to represent the rotation depends on the depth of the absorption 
profile. Smaller angular steps were not found to improve the fit 
to the optical depth r. We found that the resulting optical depth 
T is strongly affected by the gradient of the continuum's line. 
We finally look for the line at a certain rotational center whose 
corresponding linear combination of extinction profiles /fext;/ re- 
produce the optical depth t the most successfully, i.e., with the 
smallest reduced chi-square;if^ (Eq.[3]l. We use the fitting routine 
presented in Schegerer et al. (120061 ) . Representative continuum 
lines are drawn in the spectra shown in Figures |2] and [3] 

The determinations of both the continuum and the ice com- 
ponents that effectively contribute to the absorption profile were 
tested by analyzing the spectrum of the YSO GL989. Its L 
band spectral distribution was previously investigated by Dartois 
& d'Hendecourt (2001). In contrast to our measurements, the 
source was formerly observed with the ISO satellite (Infrared 
Space Observatory) where the spectrum at short wavelengths 
could be retained. But when this spectrum is cut at short wave- 
lengths up to a wavelength of 2.9 /im, the resulting spectrum 
corresponds to the NAOS-CONICA spectra that we acquired 
during different weather conditions. The results of the analy- 
sis of the 3 jjm spectra of GL 989 with and without the short- 
wavelength interval are finally compared. Potential differences 
can then be considered as an approximation of the systematic er- 
ror caused by the missing interval of the NAOS-CONICA spec- 
tra up to 2.9 yum. We note that Dartois & d'Hendecourt (1200 1 1 ) 
used a Planck function to determine the underlying continuum 
of GL 989 assuming a temperature of T = 500 K. Furthermore, 
they considered extinction profiles of ice grains with a continous 
size distribution. 

Figure|5]shows the results. Independent of the rotational cen- 
ter used, the fitting results obtained from the spectra that are 
truncated at 2.9 yum differ by < 15% from the results derived from 
spectra for which the lower end of the L band is considered. The 
mass contribution of crystallized ice can be neglected in all spec- 
tra. The mass contribution of amorphous, large grains depends 
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2JB 3.0 3.2 3.4 3.6 U 2.8 3.0 3.2 3.4 3.6 3J 




»(»«"] 



Fig. 2. Extracted and photometrically calibrated L band spectra of the source YLW 16 A for slit orientations parallel and orthogonal 
to its rotational axis. The counting rates of 4 x FWHM of the (spatial) widths of the spectra were determined. The dotted lines 
represent the "best-fit" spectral continua that we assume in our modeling approach. 



strongly on the rotational center used. Both the non-truncated 
and truncated spectra exhibit only marginal evidence of large 
water ice grains (<1%) when a rotational center at 3.8 /^m is 
used for the linear continuum. In contrast, the contribution of 
amorphous, large grains with a size of 0.5 //m is >85% when a 
rotational center at 3.6/im is used. The fitting results deviate by 
between x^ - 0.0015 and x^ - 0.003 from the measurement 
(Eq. [3]). The deviation from the measurement is larger for fits 
with the rotational center at 3.6 jum. 

5.2. Modeling results 

Table |2] and Figures |6] and |7] summarize our modeling results. 
The quality of each fit is characterized by x^ (e.g.. Press et 
al. [T986j with 



x'- 



1 



N.i-f 



Y,(fT°''\A,) - r^'HA-df 



(3) 



The data with the largest x^ have the lowest signal-to-noise ra- 
tios. In Eq. |3] A^^ is the number of data points, / the degree of 
freedom, and T^°^''^{A) and Ty'^^'^iA) are the modeled and mea- 
sured fluxes, respectively. 

The features of additional compounds (e.g., of ammonia hy- 
drate) is superimposed on fractions of the ice absorption band at 
a wavelength of ~3.4yum. The spectral region between 3.25 //m 
and 3.65 jum is not considered by our modeling approach (dotted 
curves). 



In contrast to our results for spectra with the parallel slit ori- 
entation of the slit, a larger amount of crystalline grains (28%) 
contributes to the spectrum of the orthogonal slit orientation (3'^'' 
line in the Table|2|i. We point out that this spectrum can be fitted 
by a larger amount of large, i.e., 0.5//m-sized amorphous grains, 
if crystallized grains are excluded from the fit. As the fit then 
worsens by 35%, this finding can be considered as evidence of 
crystalline dust grains, as the maximum of the extinction profile 
of larger amorphous grains shifts to longer wavelengths similar 
to the shift of the extinction profile of crystalline grains. 

In Fig. 121 the spectra 1 to 7 correspond to different regions 
within the source that could be spatially resolved by our obser- 
vations with NAOS-CONICA for the orthogonal slit orientation 
(Sect. [3]). Because of the too noisy data, we do not consider the 
optical depth derived from the observation where the spectro- 
scopic slit was orientated parallel to the rotational axis. As we 
obtain tighter fits to the derived optical depth, i.e., smaller x"^^ 
we focus only on the results obtained from fits where the con- 
tinuum's line is rotated around 3.8 jum but without forgetting the 
corresponding results obtained for a rotational center at 3.6 yum. 

The spectra 3, 4, and 5 that were derived from the most cen- 
tral region of the source around the photocenter appear to contain 
the largest contribution of non-evolved, i.e., amorphous, small 
grains, while the spectra from the outer source regions exhibit 
a greater contribution from crystallized water ice. These crys- 
talline, peripheral grains also contributes to the spectrum where 
the counting rates of 4 x FWHM were summarized (3'^'' line in 
the Table Ell. 
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Considering Fig.[3]and the flux scales used, it is conspicuous 
that different gradients are found for fitting the continua of spec- 
tra that originate in different source regions. However, if steeper 
gradients were used to fit continua of spectra originating in outer 
source regions, even a larger contribution of crystallized grains 
would be found. On the other hand, a shallower gradient for the 
continua derived in the more central source regions would favor 
a larger contribution of small amorphous grains. In any case, the 
fit is poorer when gradients other than those are used. 

The fit results obtained for spectra in which a rotational 
center at 3.6 yum is used, have qualitatively the same increas- 
ing contribution of crystalline grains towards the outer regions 
(otc - 57%, 0.25%, 98% at the regions 1, 3, and 7, respectively). 
However, as already mentioned in Sect. 15.11 the contribution of 
amorphous, large (i.e., 0.5;um-sized) grains has increased at the 
expense of amorphous, small grains, in particular in the more 
central regions 4, and 5, where on average m^^x a 88%. 

The column density A^a of the absorbing ice material can 
be determined by an integrating of the extinction profiles of 
amorphous water ice that are fitted to the optical depth t(A). 



According to Whittet (120021 ). we obtain 



A?A(H20[ice]) = r 



^■^f"^ TiA) dA 



A A^' 



(4) 



The integral extinction cross-section A of amorphous water ice at 
a temperature of r = lOKisA = 2.0xlO"'**m/molecule(Hagen 
& Tielens I1981bl l. As the O-H stretching mode is very strong 
and dominates over other species (Tielens et al. 119841 ), the 3 fiia 
band provides an accurate water ice column density. Teixeira & 
Emerson ( 1999) found the following empirical relation between 
the column density A^A(H20[ice]) and the visual extinction Ay: 



A^A(H20[ice]) = (Av - 2.1) x 10^'m 



-2 



(5) 



When considering spectra obtained from different source re- 
gions, the quantities A^a and Ay are higher by 18% and 16%, 
respectively, in spectra obtained from outer than for those from 
the inner regions. 
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Table 2. Relative mass contribution of amorphous-small (m-^y. 0.1 jum and 0.3 /im), amorphous-large (m^y. 0.5 /im, 0.8 |/m, and 
1.5 //m) and crystallized {111^: all particle radii) water ice grains derived from our modeling approach by using a rotational center 
at 3.8 ;um. The position angles (PA) of the spectroscopic slit, measured from North to East, are additionally listed. The regions 1 
to 7 represent spectra that originate in source regions with widths of 40 AU, 30 AU, 20 AU, 20 AU, 20 AU, 30 AU, and 40 AU. The 
spectra 3-5 originate closer to the photocenter than the spectra 1, 2, 6, and7. In contrast, for the first three exposures listed, the 
counting rates within 4 x FWHM of the spatial dimension are summarized. The quantities A^a and Ay represent the column density 
of water ice particles and the visual extinction derived along the line of sight (Eqs.|4] and|5]l. 
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5.3. Water ice growth in fhe disks around TTauri stars? 

The spectra of central source regions derived from a rotational 
center of the continuum's line at 3.8 /^m show profiles of amor- 
phous water ice grains that exhibit only weak hints of grain 
growth (Table |2]i. However, when considering the same results 
for a rotational center at 3.6 //m, the contribution of the ab- 



sorption profiles of 0.5;um-sized grains is significantly higher 
although the fit worsens by 24% up to 32%. From Eq. |2] we 
find that scattered radiation is present only for grain sizes of up 
~0.5fim. If we replace the extinction profiles of large grains 
(a > 0.5 yum) with the corresponding absorption profiles, the 
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spectroscopic contribution of large grains increases by only a 
few percent. 

In previous studies, water ice grains larger than 1 /im were 
not detected in YSOs. Thi et al. ( 120061 ) point to broadened wa- 
ter ice profiles in the spectra of several YSOs in the molecu- 
lar cloud Vela indicating ice grain growth. There is at least one 
among five objects in their sample where the absorption profile 
of 0.5 jum-sized, amorphous ice grains can be used for model- 
ing. Whether this result and the other evidence we have found 
indicate ice grain growth remains to be discussed in the context 
of the following issues. As shown in Schegerer et al. ( 120061) . in- 
creasing the porosity of grains with constant radii broadens the 
lOyum-silicate feature. We assume that porous water ice grains 
show similar behavior but this assumption can only be confirmed 
by future, theoretical studies. Furthermore, it is often assumed 
that ice grains have cores, such as those of 0. 1 yum-sized, amor- 
phous silicate (e.g., Jones & Merrill 1976). The core serves as 
a seed where water is adsorbed. As shown in Fig. [8] a growth 
in the silicate core also results in a broadening of the 3jum- 
absorption band (e.g.. Smith et al ll988l ). A broadening of the 
extinction profile can also be evoked by the shape of the contin- 
uum curve derived. In our study, the continuum is represented 
by a straight line, while the Planck function or spline functions 
are used in other studies. The Planck function, which is con- 
cave down, implies that there has been a broadening of the water 
ice band at long wavelengths, which automatically indicates ice 
grain growth. Finally, it is not clear to us why only 0.5 //m-sized 
grains are found when large dust grains are used in the fit. It has 
to be clarified whether the ice grain growth to sizes a > 0.5 fim 
is physically prevented in YSOs or if this finding can be ascribed 
to the fitting approach used. 

We conclude that an unambiguous identification of grain 
growth is difficult to make. In contrast, the discrimination be- 
tween crystallized and amorphous grains should be easier as ar- 
gued by Dartois & d'Hendecourt ( 120011 ). as the extinction profile 
of crystallized ice narrows. 



5.4. Crystallized water ice in the disks around TTauri stars? 

Because of the low temperatures and weak shielding of the in- 
terstellar media against hard irradiation, we assume that there is 
a negligible contribution of crystallized water ice when YSOs 
start to form from molecular clouds. The finding of crystallized 
water ice in YLW 16 A is now discussed. 

As described in Sect.|3] information about the spatial ice dis- 
tribution within the circumstellar environment can be extracted 
from the spectra measured with NAOS-CONICA. Instead of 
summarizing the intensities of the aperture, the Gaussian func- 
tion that is fitted to the intensity distribution in each row of 
the image is symmetrically divided in steps of 20 AU, 20 AU, 
30 AU, and 40 AU starting at the photocenter of the source. 
Independently of the rotational center used, we mainly find 
amorphous water ice close to the photocenter (ma > 88%). 
However, for increasing distance from the photocenter, the mass 
contribution of crystallized ice strongly increases (Table |2] 
Figl3- The sketch in Fig.[TO]shows a circumstellar disk where 
the location of amorphous and crystallized water ice grains are 
plotted considering our fitting results. We assume that the crys- 
tallized water ice component has its origin in the upper disk lay- 
ers at large disk radii, where water ice does not evaporate and is 
protected by the surrounding material of the circumstellar enve- 
lope from iiTadiation. Because of the higher temperatures, water 
ice at smaller radii close to the snowline is assumed to be crys- 
tallized inside the more embeded disk layers. The observation 
of the outer parts of the optically thick disk also explains the 
higher column density A^a inferred from spectra that originates 
in the outer regions (Table |2|). Again, this finding does not de- 
pend on the rotational center used. The photocenter of the spec- 
tra is assumed to represent the optically thin bipolar envelope 
of the source observed with NICMOS/HST (Allen et al. 2002] ). 
Because of the low effective shielding in this central region, crys- 
tallized water ice is assumed to be destroyed by irradiation. 

Figure |9] shows the normalized optical depths derived from 
the depth of the peripheral region 7 {gray curve) that is overlaid 
on the mean of the depths of more central regions. The optical 
depths of the more central regions are lower than those of the 
peripheral regions. We note that if only amorphous grains are 
used in our deconvolution approach, the fitting worsens by up to 
14% according to the parameter ;\f^ (Eq.[3j Table|2]) derived. 

The YSOs that were additionally observed with NAOS- 
CONICA and whose observations are presented in the Appendix 
of this publication, do not show any evidence of crystallization. 
However, the spectra of these objects contain the water ice ab- 
sorption feature. The question of whether these features orig- 
inate in the circumstellar environment of these objects cannot 
unequivocally be answered here. 



6. Summary 

We have presented an L band spectroscopic observation per- 
formed with NAOS-CONICA of a YSO that is strongly inclined. 
These systems allow the study of optical thick circumstellar re- 
gions that shield water ice from hard irradiation and evaporation. 
The absorption feature of water ice at 3.08 //m was identified in 
the spectra with optical depths of between 1.8 and 2.5, depend- 
ing on the region observed. The optical depth, the water-ice col- 
umn density, and the visual extinction are higher in the more pe- 
ripheral regions at greater distance from the photocenter of our 
observations. The optical depths were derived from the spectra 
by assuming that a straight line represents the continuum flux. 
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Fig. 9. Optical depths derived from the central (black curve) and 
peripheral (gray curve) regions of the YSO observed. The opti- 
cal depths of the region 3, 4, and 5 were averaged to obtain the 
black curve while the depth from the region 7 was used to ob- 
tain the gray curve. For comparison, the curves are normalized 
at the optical depth at -2.9 yum. A larger amount of crystalline 
ice grains contribute to the narrow gray curve with a maximum 
at 3.075 ;um. Almost only amorphous dust grains contribute to 
the black curve derived from the more central region. 



Considering the extinction profiles of amorphous and crys- 
tallized water ice grains with sizes ranging from 0.1 yum up to 
1 .5 fiui, the derived optical depths were deconvolved using a fit- 
ting routine that was already presented in Schegerer et al. (I2006I I. 
The optical depths of the spatially unresolved spectra were found 
to be dominated by the extinction profile of small, amorphous 
(i.e., non-evolved) water ice grains. However, in spectra for 
which the T Tauri object YLW 16 A could be spatially resolved, 
crystallized water ice could be found in possible outer disk re- 
gions, i.e., in disk layers at radii >30 AU in projection. The pu- 
tative bipolar envelope of this source was not found to exhibit 
any evidence of crystallized material (Fig.fTOll. A growth in the 
size of the water ice grains cannot unambiguously be determined 
because the corresponding modification of the water ice feature 
may also be either caused by an increase in the porosity and/or a 
growth of the particle core. 
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Fig. 10. Sketch for a strongly inclined circumstellar disk such as YLW 16 A. We assume that amorphous, non-evolved water ice is 
located in the circumstellar disk and in the bipolar envelope. The assumed location of crystallized water ice in the disk is close to 
upper disk layers, at disk radii that are larger than that of the snowline. The observer and the spectroscopic slit are also sketched. 
The orthogonal slit orientation is shown. Radiation from amorphous and crystallized ice grains are represented by dashed and solid 
arrows, respectively. 
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Appendix A: Additional measurements 



Table lA. ll provides details of additional observations with NAOS-CONICA. The quantity Ttot is the total exposure time. The standard 
stars HR 6070 and BS 7330 were observed directly before or after the observation of each target. The brightness of the standard star 
BS 7330 is known from the standard star catalog of Leggett et al. C2003J . The observation of the standard star BS 7330 helped to 
remove the telluric absorption features from the spectrum of VV CrA and its subsequent photometric calibration, while the object 
HR 6070 was used as a standard star for all remaining objects. 



The L band brightness of the objects CRBR 2422.8-3423 and VVCrA A are found in the Gezari-catalog ( fT999] l. The L band 
brightness of the infrared companion VV CrA B is not known. The FWHM of the projected diameter of CRBR 2422.8-3423 could 
be estimated using a K band image (Pontoppidan et al. '20051 The FWHM of the diameter of VVCrA A and VVCrAB were 
derived using interferometric measurements in N band (Przygodda 2004). The flux and the diameter of the object Elia 2-23 in the L 
band were taken from Haisch et al. (2002). The NICMOS/HST-image of GY92244 in the H band (Aflen et al. l2002t aUows us to 
approximate the diameter of the corresponding disks in L band. The L band brightness of this object was estimated by considering 
known H band magnitudes and (H-K) and (K-N) colors of additional YSOs in the same star-forming region (Allen et al. l2002l) . 



Analogous to YLW 16 A, all targets were observed twice, with two orthogonal orientations of the slit: parallel (p) and orthogonal 
(o) to a predefined direction. If the position angle of the disk was not known, the orientation of possible companions were considered 
in the subsequent selection of the slit orientation. A slit length of 28" allowed us to observe the companions of VVCrA A (i.e., 
VVCrAB at an angular distance of 1.69"; Przygodda '2004^ and Elia 2-23 (Elia 2-21 at an angular distance of 10"; Haisch et 
al. |2002J in the parallel orientation of the slit, simultaneously. The data reduction of these observations is described in Sect.[3j 



Fi gures [ATI and [Am show the resulting L band spectra of these targets. Considering the derived flux, some sources exhibit an 
intrinsic variability in L band (Leinert et al. 1200 11 1. For instance, it is known that the NIR brightness of VV CrA "conspicuously" 
varies within a period of few years (Graham 11 9921 Chen & Graham[T993j Koresko et al. 119971 Przvgodda l2004!) . 



The result of our observations of VV CrA AO conflicts with a spectroscopic measurement performed on data acquired using the 
4m-telescope of the CTIO (Cerro Tololo Inter- American Observatory; Chen & Graham [T993]l. The absorption band of VV CrA A 
detected by our observation is weak, while the observation of Chen & Graham (1993J does not show any absorption ice band at 
all. This difference may be attributed to the factor-of-two higher spectral resolution power of NAOS-CONICA and its higher spatial 
resolution power by a factor of ~ 1 1 . Therefore, observations with NAOS-CONICA allow us to observe more central regions that are 
strongly influenced by both the star and disk. The water ice absorption band of the infrared companion is more pronounced. This 
companion, that is strongly variable in the infrared wavelength range (Koresko et al. 1997\ exhibits a broad absorption band. Our 
analysis of the observations of this source agree with previous measurements with the 4m-telescope at the CTIO. 



^ The main component of the binary in N band is the south-west component in the binary system. The north-east component is brighter in L 
band. 
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Table A.l. Overview of additional observations with NAOS-CONICA. Coordinates, L band magnitude, FWHM, date, observing 
time (UT), and the mean airmass (AM) during the observations are listed. 



object 



Q'(J2000) (5(J2000) L [mag] FWHM ["] 



date 



UT 



AM 



April 2"", 
Apr 2'"', 

April 2'"', 
April 2"'', 

June 22'"', 
June 22'"', 



HR6070 


161818 


-28 3646 


4.807 


Elia2-23 (p) 


16 2624 


-2424 50 


5.9 


HR6070 


161818 


-28 3646 


4.807 


Elia2-21 


162624 


-2424 50 


8.0 


Elia2-23 (o) 


1626 24 


-2424 50 


5.9 


HR6070 


161818 


-28 3646 


4.807 


GY92 244(p) 


1627 17 


-2428 58 


8.8 


HR6070 


161818 


-28 3646 


4.807 


GY92 244(o) 


16 27 17 


-2428 58 


8.8 


HR6070 


161818 


-28 3646 


4.807 


CRBR 2422.8-3423 (p) 


16 27 05 


-24 36 31 


9.7 


HR6070 


161818 


-28 3646 


4.807 


CRBR 2422.8-3423 (o) 


16 27 05 


-2436 31 


9.7 


HR6070 


161818 


-28 3646 


4.807 


VVCrAA(p) 


19 03 07 


-37 12 50 


3.7 


BS7330 


19 2129 


-3458 57 


4.9 


VVCrAA(o) 


19 03 07 


-3712 50 


3.7 


BS7330 


19 2129 


-3458 57 


4.9 


VVCrAB 


19 03 07 


-37 12 50 


6.2 


BS7330 


19 2129 


-3458 57 


4.9 



1 



2006 
2006 

2006 
2006 

2006 
2006 



02-7: 
45-9: 
02-7 : 
45-9: 

00-0: 
45-0: 



04 

34 

04 

34 

30 

47 



May 27"', 2006 2:48-4:00 

May 27"', 2006 4:24-4:26 

May 27"', 2006 5:25-6:02 

May 27"', 2006 6:22-6:24 



1 August 12"', 2006 2 : 01 - 2 : 57 

Aug 12"', 2006 3 : 15-3 : 17 

1 August 19"', 2006 : 14 - : 45 

- August 19"', 2006 1 : 11 



1 : 13 



July 27'*, 2006 6:28-6:30 

July 4"', 2006 5:17-6:03 

July 4"', 2006 6 : 20 - 6 : 22 

July 27"', 2006 5:23-6:01 

July 27"', 2006 6 : 28 - 6 : 30 



1.0 
1.0 

1.0 
1.0 

1.0 
1.0 

1.1 
1.0 

1.0 
1.1 



July 27"', 2006 5:23-6:01 1.2 



1.3 

1.0 
1.1 

1.2 



2.0 s 
72x6x2.0s 

2.0 s 
72x6x2.0s 

76x6x2.0s 
2.0 s 

54 X 18 X 2.0 s 
2.0 s 

36 X 18 X 2.0s 
2.0 s 



1.3 18 X 40 X 4.0 s 

1.5 2.5 s 

1.0 12 X 60 X 3.0 s 

1.1 2.5 s 



60x8 X 1.5 s-f 

-1-16 X 5 X 5.0s 

2.0 s 

60x 12x 1.0 s 
2.0 s 



60x8x 1.5 s-i- 
-1-16 X 5 X 5.0s 
1.3 2.0 s 



Table A.2. Column density A^a of water ice along the line of sight to the observer (Eq.|4]i. 



exposure 


PAT] 


A'a [10^- m-^] 


Av [mag] 


ra.;s [%] 


ra.;i [%] 


me [%] 


X 


Elia2-23(p) 


53 


0.52 


7 


- 


- 


- 


~ 


Elia2-23(o) 


45 


0.26 


5 


100 


- 


- 


0.04 


Elia2-21 (p) 


53 


0.66 


9 


- 


- 


- 


~ 


GY92244(p) 


178 


1.8 


20 


>99 


< 1 


< 1 


0.08 


GY92244(o) 


88 


2.1 


23 


>99 


- 


< 1 


0.1 


CRBR 2422.8-3423 (p) 


53 


3 


35.1 


- 


- 


- 


- 


CRBR 2422.8-3423(0) 


45 


3 


29.1 


91 


< 1 


8 


0.2 


VVCrAB(p) 


134 


0.66 


9 


>99 


- 


< 1 


0.05 


VVCrAA(p) 


134 


0.13 


3 


- 


- 


- 


- 


VVCrAA(o) 


43 


0.030 


2 


- 


- 


- 


- 



Owing to the noisy data from which only parts of the telluric features could be removed, processed water ice, i.e., crystallized 
and/or grown ice grains cannot unambiguously be found. Residual telluric features also hamper the determination of the column 
densities A^A(H20[ice]). Figure IA3] shows the optical depths of the objects for which the determination of the column density and 
a subsequent modelling were possible, given the residual telluric features. Table lA.2] lists the coiTesponding results. 

The detection of crystallized water ice in the spectra of CRBR 2422.8-3423 cannot definitely be confirmed. Studying the derived 
profiles in Smith et al. (1989) and Dartois & d'Hendecourt (2001), a signal-to-noise ratio greater than ~4 is required to prove crys- 
tallized water ice in the spectra, i.e., to distinguish the extinction profiles of crystallized and amorphous ice grains. Our measurement 
of CRBR 2422.8-3423 is obtained from data with a signal-to-noise ratio of ~6. However, the profile of CRBR 2422.8-3423 cannot 
be used because its absorption is saturated (Fig. IA.lt . 
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Fig. A.l. Extracted and photometrically calibrated L band spectra of further targets. 
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Fig. A.2. Continuation of Fig. lA.l 
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Fig. A.3. Optical depth t(X) that is derived from the spectra of Fig. [ATI and Fig. lA.2l The optical depth of the spectra Elias 2-23 (p), 
Ehas2-23 (o), Elias2-21, CRBR 2422.8-2423 (p), VVCrA (p), and VVCrA (o) are not modeled because of low signal-to-noise 
ratios and remaining telluric features that could not be sufficiently eliminated in data reduction. 



